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We calculate in perturbative QCD the transverse momentum broadening of vector boson pro- 
duction in high energy nuclear collisions. We evaluate the effect of initial-state parton multiple 
scattering for the production of the Drell-Yan virtual photon and W/Z bosons. We calculate both 
the initial- and final-state multiple scattering effect for the production of heavy quarkonia and their 
transverse momentum broadening in both NRQCD and Color Evaporation model of quarkonium 
formation. We find that J/tp and T broadening in hadron- nucleus collision is close to 2Ca/Cf 
times the corresponding Drell-Yan broadening, which gives a good description of existing Fermilab 
data. Our calculations are also consistent with RHIC data on J/i/> broadening in relativistic heavy 
ion collisions. We predict the transverse momentum broadening of vector boson (J/tp, T, and W/Z) 
production in relativistic heavy ion collisions at the LHC, and discuss the role of the vector boson 
broadening in diagnosing medium properties. 
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I. INTRODUCTION 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory has been providing experimenters 
with colliding beams of heavy nuclei at relativistic energies as high as 100 GeV per nucleon. Early RHIC data 
from central gold-gold collisions strongly indicate the creation of dense and hot QCD matter with very unusual and 
interesting properties [1[. Characterized by its opacity to jets and other evidences, the new QCD matter shows an 
unprecedented energy density well above the critical value predicted by lattice QCD for establishing quark-gluon 
plasma (QGP), a weakly coupled gas of quarks and gluons 0. On the other hand, the data from the observation 
of an unexpectedly large flow indicate that the hot QCD matter interacts strongly with itself and behaves like an 
almost ideal liquid with low shear viscosity rather than a gas [l|. In this paper, we investigate the role of transverse 
momentum broadening of heavy vector boson production in identifying the formation of new QCD matter and in 
diagnosing medium properties. 

Hadronic production of a heavy vector boson, such as a virtual photon with a large invariant mass, a heavy 
quarkonium, and a Z (or a W) boson requires a short-distance hard partonic collision with a large momentum 
transfer because of the large vector boson mass, and is dominated by the scattering of two partons from two incoming 
hadron beams (or from a beam and a target). The short-distance creation of the vector boson or the heavy quark 
pair could be perturbatively calculated in Quantum Chromodynamics (QCD), while the hadronization of the quark 
pair to a physical quarkonium offers unique perspective into the formation of QCD bound states. Heavy quarkonium 
production in high energy nuclear collisions is of special interest in diagnosing the properties of nuclear medium 
because of the two-scale nature of the creation process. The creation of the heavy quark pair takes place at such 
a short time that it is unlikely to interfere with the dynamics of nuclear medium which is effectively frozen. On 
the other hand, the hadronization from the heavy quark pair to a bound quarkonium could be very sensitive to the 
medium properties. A change in production rate and momentum spectrum for heavy quarkonium production from a 
proton-proton collision to a nuclear collision signals a change of characteristics of the nuclear matter and provides the 
opportunities to probe the matter's properties [H, H[. The challenge is to quantitatively extract interesting properties 
of the nuclear matter from the change. 

In high energy hadron-nucleus and nucleus-nucleus collisions, the interaction between the produced heavy quark 
pair and the nuclear medium (referred to as a final-state effect), as well as possible scattering between the incoming 
colliding partons and the nuclear matter before the hard collision to produce the pair (an initial-state effect), can both 
change the heavy quarkonium's production rate and its momentum spectrum. The effect of the final-state interaction 
depends on the hadronization mechanism - how a produced heavy quark pair becomes a bound quarkonium. A color 
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singlet heavy quark pair with a small color dipole moment is less likely to interact with the nuclear medium than 
a pair in a color octet configuration [f| 0]. That is, the final-state effect could be an ideal probe for exploring the 
nonperturbative formation mechanism in heavy quarkonium production. But, the observed nuclear dependence is 
sensitive to not only the final-state effect but also the initial-state effect as well as the quantum interference between 
them. Due to the local hard collision that produces the heavy quark pair, the quantum interference between the final- 
state and initial-state partonic rescattering could be suppressed by the hard collision scale. To understand the effect 
of the initial-state interaction, we also study nuclear dependence in the production of the Drell-Yan virtual photon 
as well as Z (or W) vector boson. These vector bosons, if reconstructed from their leptonic decays, do not interact 
strongly once produced at the short-distance. The initial-state effect itself can be a good probe for the dynamics of 
partonic rescattering when a fast parton passes through the nuclear matter. 

Partonic multiple rescattering in a nuclear medium before as well as after the hard collision can modify the distri- 
bution of the vector boson's transverse momentum qT- Each rescattering is likely to change the momentum spectrum 
by an order of the typical transverse momentum of the partons inside the nuclear matter, which is significantly softer 
than the momentum exchange in the hard collision. Therefore, the change to the transverse momentum spectrum, 
da/dqj>, should be most significant when qr is relatively small. But, the rescattering effect on the low qx spectrum 
is unlikely to be calculable in perturbative QCD @, On the other hand, an averaged transverse momentum 

square of the produced heavy vector boson, 



fa) = I ddrqr^-/ J dqrp dq% ' (1) 



is much more inclusive and perturbatively calculable if the qr is integrated over a wide range 
accumulative change from the rescattering to the averaged transverse momentum square - transverse momentum 
broadening, A(<?|,) = {qj)\AB — (Sr)U/u defined as a difference between the calculable transverse momentum square 
in nuclear collision and that in hadron-hadron collision, should be calculable too 0, @ • 

The Drell-Yan transverse momentum broadening A(<7j,)dy has been studied in perturbative QCD by evaluating 
coherent partonic rescattering diagrams between incoming (anti)quark and the cold nuclear matter [13, El, EH- The 
rescattering leads to the broadening of the Drell-Yan dilepton's transverse momentum distribution. The broadening 
was shown to be proportional to the target size or to have the A 1 / 3 -type nuclear dependence. The calculated nuclear 
dependence was found to be consistent with both Fermilab and CERN data [HI]. On the other hand, as shown in 
Ref. [l5j . there have been difficulties in understanding the same broadening in the production of heavy quarkonia 
(3/tp and T). At the leading order of perturbative calculation, the Drell-Yan dilepton production is dominated by a 
quark- antiquark annihilation subproccss, while the heavy quarkonium production in hadronic collision is dominated 
by a gluon-gluon fusion subprocess. If one neglected final-state interaction in both Drell-Yan and heavy quarkonium 
production, one would expect that the ratio of broadening of heavy quarkonium over Drell-Yan is close to the ratio of 
the multiple scattering effect of a gluon over that of a quark (or an antiquark) inside a nuclear medium [l^ |. At the 
lowest order, the ratio is approximately equal to the ratio of color factors of the lowest order gluon and (anti)quark 
rescattering, C A /C F = 9/4, with C A = N c = 3 and C F = (N^ - l)/(2iV c ) = 4/3. Although the data on heavy 
quarkonium broadening in hadron-nucleus collisions shows the expected A x l -t ype nuclear dependence, the ratio to 
the Drell-Yan broadening could be as large as 5, twice of the naive expectation [15| . 

Recently, Johnson et al. published in Ref. [17j a new analysis of Drell-Yan broadening, A(q^)oY, based on the 
ratio of Drell-Yan transverse momentum distribution on heavy nuclear targets over deuterium target, Ra/d- By 
using the chain relation, Ra/d — Ra/Bc Rbc/Di the new data from Fermilab E866 experiment on the ratio Ra/Be 
with A = Fe,W over beryllium (Be), and a theoretical formula for Rse/Di Johnson et al. showed the result of 
their analysis in Fig. 5 of Ref. [13] and concluded that the observed broadening is about twice as large as the one 
published previously [l5j . which implies that the heavy quarkonium broadening is just a factor of Ca/C f larger 
than the Drell-Yan broadening and is consistent with the naive expectation. The factor of two enhancement of the 
Drell-Yan broadening represents a difference of more than three standard deviations between the central values of two 
different analyses in terms of the error size of previous analysis or a difference of two standard deviations if one uses 
the much larger error of the new analysis. It will be interesting to see if this new analysis changes the broadening of 
J /tp and T as well. If the ratio of the heavy quarkonium broadening over Drell-Yan broadening is indeed consistent 
with the naive expectation, we would have to conclude based on the discussion on final-state interaction in this paper 
that the heavy quark pair in heavy quarkonium production is mainly produced in a color singlet state, and the color 
of an octet pair will have to be neutralized at a very short distance. This conclusion would have a very important 
impact on the physics of J ftp suppression in nuclear medium. 

Since a heavy quarkonium is unlikely to be formed at the same time when the heavy quark pair was produced [l8| , the 
final-state interaction between the heavy quark pair and the nuclear medium could generate additional broadening. 
Since the broadening is evaluated from the transverse momentum distribution of observed heavy quarkonia and 
normalized by the total production rate, as defined in Eq. (TTJ), we will not discuss the impact of final-state interaction 
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on the suppression of heavy quarkonium production in this paper [4, 5]. Since the final-state interaction is sensitive 
to the non-perturbative formation mechanism of the bound state, we calculate the broadening of heavy quarkonium 
production in both non-relativistic QCD (NRQCD) model and Color Evaporation model [2(|. The often used 
Color Singlet model [2l[ can be interpreted as a special case of the NRQCD model. In the NRQCD model, the non- 
perturbative dynamics for a heavy quark pair to form a bound quarkonium is organized through matrix elements of 
operators that are characterized by an expansion in the relative velocity of the pair and the pair's rotational and color 
quantum numbers. On the other hand, in the Color Evaporation model, all heavy quark pairs with invariant mass 
less than the mass threshold of producing a pair of open flavor heavy mesons have the same probability to become 
a bound quarkonium regardless the pair's rotational and color quantum numbers. Rescattering of the heavy quark 
pair in nuclear medium could change the pair's rotational and color quantum numbers. Therefore, these two models 
could lead to different predictions for the nuclear dependence of heavy quarkonium production in nuclear collisions. If 
the difference is significant, an accurate measurement of the nuclear dependence could provide important information 
on the hadronization. We find that the heavy quarkonium broadening calculated in these two models have different 
analytical expressions. But, numerically, these two models predict a very similar result for the broadening of J/ip and 
T in hadron- nucleus collision. The calculated broadening is close to 2Ca/Cf, which is consistent with Fermilab data 

MM- 

We also investigate the nuclear dependence of the averaged transverse momentum square in relativistic heavy ion 
collisions. If all soft gluons of heavy ion beams are stopped to form the hot dense medium, the final-state interaction 
between the slowly expanding or almost stationary medium and the fast moving heavy quark pair (or a quarkonium) 
of transverse momentum qr is unlikely to broaden the qr spectrum. Instead, the final-state interaction is likely to 
slow down the heavy quarks due to the energy loss [2j|, and could reduce the averaged transverse momentum \qj) 
[24| . On the other hand, the initial-state interaction should behave very similar to that in hadron- nucleus collision 
up to a possible difference in the medium density. In order to separate the effect of final-state interaction from that 
of initial-state interaction, and independently study the initial-state interaction and extract the medium density, we 
calculate the broadening of Z (as well as W) bosons in relativistic heavy ion collisions at the Large Hadron Collider 
(LHC). Since there is effectively no final-state interaction for the production of Z (or W) bosons extracted from their 
leptonic decay, the broadening is an ideal probe of the density of nuclear medium in the early stage of the collision [25| . 
If one would be able to reconstruct Z (or W) bosons from their hadronic decay, the final-state interaction between 
the decaying jets and the medium could lead to an apparent mass shift for the vector boson and provide additional 
tools to extract the medium's properties. 

The rest of our paper is organized as follows. In Sec. |H] we review the perturbative QCD calculation for the 
transverse momentum broadening in Drell-Yan production to set up the notation and terminology. In Sec. IIII1 
we derive the transverse momentum broadening of heavy quarkonium production in hadron-nucleus collision. We 
calculate both initial-state and final-state multiple scattering. We evaluate the transverse momentum broadening in 
both NRQCD model and Color Evaporation model. We also discuss the broadening of heavy quarkonium production 
in nucleus- nucleus collisions. In Sec. IIV1 we calculate the broadening of Z (as well as W) boson production in both 
hadron-nucleus and nucleus- nucleus collisions at the LHC. We argue that the transverse momentum broadening of 
Z (or W) bosons that are reconstructed from their leptonic decay channels is an excellent probe for initial-state 
rescattering and the medium density at an early stage of relativistic heavy ion collisions. We present our numerical 
study of the transverse momentum broadening of vector boson production in Sec. El We discuss the extrapolation of 
the non-perturbative matrix elements from the fixed target energies to collider energies. We compare our calculations 
with data from both fixed-target experiments at Fermilab, and collider experiments at RHIC. We also predict the 
broadening at the LHC energy. Finally, in Sec. IVI1 we conclude with a summary that suggests directions for future 
work. 



II. TRANSVERSE MOMENTUM BROADENING IN DRELL-YAN PRODUCTION 

The Drell-Yan production of a massive pair of leptons in hadronic collision has been a potent probe of short- 
distance dynamics in strong and electroweak interactions. The lack of interaction between the produced leptons and 
the hadronic matter makes the Drell-Yan massive dilepton production in hadron-nucleus and nucleus-nucleus collisions 
an ideal hard probe of initial-state partonic scattering in nuclear matter. At the fixed target energies, the Drell-Yan 
process is dominated by the production of a virtual photon of invariant mass Q that decays into the measured lepton 
pairs, while the production is dominated by the Z boson at the LHC energies. In nuclear collisions, it is very likely 
that the energetic incoming parton can have several scatterings with soft partons inside the nuclear matter before the 
hard collision to produce the vector boson. Such initial-state multiple scattering could induce more soft radiation from 
the incoming parton and broaden the parton's transverse momentum. The additional parton transverse momentum 
at the hard collision leads to the broadening of the observed vector bosons. 
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The transverse momentum distribution of the Drell-Yan vector boson production at high qx is calculable in per- 
turbative QCD [2(1 But, the distribution at low -C Q, which is sensitive to the soft rescattering, is not 

perturbatively calculable unless the resummation of large Sudakov logarithms dominates the shape of the low qr 
spectrum [ToL l28l . [2^ . [30j | . Each soft rescattering can only change the incoming parton's transverse momentum by the 
amount close to typical virtuality of partons inside the nuclear matter, and is too soft to be perturbatively calculable. 
Furthermore, the potential interference between the rescattering and the parton shower, which is responsible for the 
large contribution of Sudakov logarithms at low qx, could complicate the resummation of the logarithms and lead 
to even less control on the low qx spectrum. On the other hand, the averaged transverse momentum square, (q T ) 
defined in Eq. |T]), is much more inclusive. If we integrate over all kinematically allowed qx, the (q T ) depends on 
only one single hard scale, Q, the mass of the vector boson and is perturbatively calculable. The large logarithmic 
contribution to the ^-distribution from the power of ln(Q 2 /q T ) is suppressed by the q T weight. 

The transverse momentum broadening, A(q T ) = (q T )\AB ~~ (Qt) \hh, which sums over the accumulative effect of many 
soft rescattering, is expressed in terms of the difference of two inclusive and perturbatively calculable quantities, and 
is therefore calculable in perturbative QCD 0, @|. The broadening of the Drell-Yan production of lepton pairs in 
hadron- nucleus collisions was first studied in terms of a non-relativisitic QED model in Ref. [l2j • It was shown that 
initial-state interactions lead to an increase in the average of the Drell-Yan dilepton's transverse momentum square 
and the increase is proportional to the length of the nuclear target. The Drell-Yan transverse momentum broade ning 
was also systematically studied in terms of perturbative QCD collinear factorization approach in a covariant gauge [131 ] 
and was further studied in Ref. [14j in a light-cone gauge. Since we calculate the transverse momentum broadening of 
the heavy quarkonium production in a covariant gauge in this paper, we briefly review the perturbative QCD collinear 
factorization approach and the covariant gauge derivation of the Drell-Yan broadening in the rest of this section. 

The cross section for the Drell-Yan process in hadron-nucleus collisions, h(p') + A(p) — > 7* («?)[—» l + l~] + X, where 
q,p' ,p are the four momentum of the virtual photon, the incoming hadron, and the nucleus (per nucleon) with atomic 
weight A, respectively, can be expanded in terms of contributions with different number of rescattering, 

OKA = °hA + a hA + ■■■ ( 2 ) 

with superscript S for single scattering, D for double scattering, and etc. A single hard scattering is localized in 
space and time, and is unlikely to provide the target length (or the ^4 1//3 -type nuclear size) enhancement to the cross 
section, although it can get a weaker nuclear dependence to the cross section from nuclear parton distributions 0]. 
The leading contribution to the broadening of the dilepton's transverse momentum square comes from the double 
scattering [131 ] . 

with the inclusive Drell-Yan cross section given by 

da hA da% A A ^ f , 1 , n f , , , s da, 



dQ 2 dQ 2 



A^Jdx' fa /h (x>) J dx<t> g/A (x) ^ , (4) 



where A is the atomic weight of the nucleus, Y] q runs over all quark and antiquark flavors, (j)q/h and (pq/A represent the 
hadron and nuclear partonic distribution functions, respectively, and da q qjdQ 2 is the lowest partonic qq annihilation 
cross section to a lepton pair of invariant mass Q. In Eq. ([4]) and the rest of this paper, we suppress all dependence on 
the factorization and renormalization scales. In Fig.[TJ we sketch the leading order Feynman diagram that contributes 
to the double scattering cross section, dcr® A . As shown in Fig. [TJ an antiquark of momentum x'p 1 from the incoming 
hadron scatters off a gluon from the nucleus (indicated by the bottom blob) before it annihilates with a quark from 
the nucleus to form a vector boson of large invariant mass, Q, which then decays into a lepton pair. The interference 
diagrams, that have both gluons in the same side of the final-state cut (the dashed line), do not contribute to the 
broadening in a covariant gauge calculation [l3T |. while they are very important in the light-cone gauge calculation 
[Til ]. It is clear from the diagram that the momentum of the observed vector boson is only sensitive to the total 
momentum from the nucleus, which is equal to a sum of the gluon and quark momentum. Therefore, the gluon 
(or quark) momentum in the scattering amplitude (the left of the dashed line) is not necessary to be equal to the 
gluon (or quark) momentum on the right of the final-state cut. This is a consequence of the fact that there could be 
an arbitrary momentum flow from the nucleus through the quark line, the internal antiquark line, and back to the 
nucleus from the gluon line without changing both initial- and final-state. To drive the double scattering contribution 
to the cross section, we need to integrate over this loop momentum for both the amplitude and complex conjugate 
of the amplitude, or equivalently, the momentum flows through those two gluons in Fig. [TJ The internal antiquark 
propagator following the gluon rescattering can be very large if the gluon momentum is very soft, and it can actually 
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diverge if the gluon momentum vanishes. But, it is easy to verify that the singularity of the internal antiquark 
propagator when gluon momentum vanishes is not pinched. The integration of the gluon momentum can be deformed 
far away from the on-shell singularity into a perturbative off-shell region at the order of the hard scale Q; and the net 
result from the integration is given by the residue of the pole of the antiquark propagator 0, [8[ . 




FIG. 1: Lowest order double scattering Feynman diagram that contributes to the broadening of Drell-Yan transverse momentum 
distribution, which shows an antiquark of momentum x'p' of incoming hadron scatters off a gluon of a nucleus (the bottom 
blob) before it annihilates a quark to produce a vector boson. 



Following the derivation in Ref . , the contribution from the double scattering diagram in Fig.[T]to the g^-moment 
of Drell-Yan cross section can be expressed as 



dlr Qt ^Q^dq 2 = ^ J dq ^ q T J dx'&q/hix') J dx dxi dx 2 d 2 k T T Aq {x, x%, x 2 , kr, 



P) 



x H(x, xi,x 2 , kr,p, q, x'p') 5(qj. - fcf.) , 



(5) 



where the matrix element T ' A q is given by the bottom blob in Fig. [TJ which includes the propagators of all quarks and 
gluons connecting to the blob, 



T- , i , I dy dy x dy 2 f d 2 y T 

T A q(x,Xl,X2,k T ,p) = I — — —J — 

1 



e ixip + y 1 e i{x-x t )p + y e -i(x-x 2 )p + y 2 e ~ik T -y T 



x 1 ( PA I A+ (j/ 2 - , T )^(0) 7 + i>(Vi )A+ (y- , y T ) \pa) 



(6) 



where the subscript "Aq" indicates that the matrix element is made of the gluon and quark field operators. The 
partonic part H in Eq. ([5]) is given by the top partonic part of the diagram in Fig. Q] with two antiquark lines traced 
with (7 -p')/2, two quark lines from the nucleus traced with (7 •p)/2, and the Lorentz indices of two gluon lines from 
the nucleus contracted by p a p 13 [13| . 

The separation of the partonic part H from the hadronic matrix element T Aq in Eq. ([5]) is not yet a gauge 
invariant factorization. The matrix element T Aq in Eq. ([6|) has an explicit dependence on the gluon field operator 
A + , which is not gauge covariant. That is, the matrix element T Aq cannot be made gauge invariant by the insertion 
of ordered gauge links between the field operators [3l|. To achieve a gauge invariant factorization, it is necessary 
to convert the gluon field operator A + in the matrix element into corresponding gluon field strength, F +a , with 
a transversely polarized Lorentz index a. We can achieve this conversion in a covariant gauge as follows. We 
first expand the kx in the partonic part, H, around fcf, = 0, because k\ -C Q 2 , and keep the first nonvanishing 
term, H(x,x\,X2,p,q 1 x'p') — H(x,xi,X2,kr = 0,p,q,x'p'), which is given by the diagram in Fig. [2] We then 
write the transverse momentum square, q\ in Eq. §5§ as k\ by taking advantage of the S(q 2 ^ — and convert 
k\ A + (y2 ,0r)A + (y~ ,yr) to F a + {y 2 -,0T)F +a (y~ ,yr) by a partial integration [l3l| . Finally, upto the power corrections 
in (fc|n)/Q 2 , we can rewrite the ^-moment in Eq. ([5]) as 



j 2 2 
dq T q T 



dQ 2 dq 2 



^ J dx' (j)q/h{x') J dxdxi dx 2 T^(x,xi,x 2 ,p) H(x,x 1 ,x 2 ,p, q, x'p) 



(7) 




(x-xi)p 




(X-Xi)p 



FIG. 2: Lowest-order double scattering diagram that leads to the factorized partonic part, H, in Eq. ((9j). 
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where Tp q is a twist-4 parton correlation function defined as 



tP(x,Xi,X2,P) = I ^— UVl i^2_ e " lP+yi e i ^ x ~ Xl ^ p+v ~ e -i(.x-w)p + v~2 



dy dy 1 dy 2 
2?r 2tt 2tt 



x \ (pa I F a + (yj )^(0) 7 + ^(yr )^ +Q (iT ) M 



(8) 



with the superscript "(/)" indicates the matrix element corresponding to the initial-state rescattering [13j . The leading 
order contribution to the partonic hard part from the diagram in Fig. [2] is fl3l ] 



H(x,x 1 ,x 2 ,p,q,x'p' 



N? - 1 



1 



1 



2tt x\ — x — if. X2 — x + ie 



d(T q q 

~d& ' 



(9) 



where two unpinched poles are from the two antiquark propagators in Fig. O Substituting Eq. to Eq. ([7]), and 
integrating over x\,X2, by taking the residues of the unpinched poles (the leading pole approximation [7J), we obtain 



/ 



2 2 MCT /iA 



dQ 2 dq 1 \ 



J2 / dx' cj> q/h (x') / dxT^x) 



d<7gg / 87T 2 a s 



dQ 2 V^c-l 



C, 



(10) 



with the measurable twist-4 quark-gluon correlation function 0, fl3| , 



2tt 



dyi dy 2 
2ir 



■0(y- -yi)0(-V2) 



x g (PA|F« + (j/ 2 -)^(0)7 + ^(2/-)^ +a (yr)bA> , 



(11) 



where the superscript "(/)" again indicates the initial-state rescattering. From Eq. ([3]), we obtain the leading double 
scattering contribution to the Drell-Yan broadening (l3j . 



T/DY 



</<T„ 



^Y^ql dx 1 <j) q/h {x') j dx(j) q/A {x) dQ , 
By using the model proposed for the twist-4 parton correlation functions 

urn 

T^(:r)=A 2 ^/% M (z), 
we can express the Drell-Yan broadening in a much simpler form [l3j 

8tt 2 q; 



T/DY 



TV 2 - 1 



A 2 A l/3 



(12) 



(13) 



(14) 



with an unknown non-perturbative parameter A 2 defined in Eq. (Ti"3"]) . The leading contribution to the Drell-Yan 
broadening in Eq. (Ti"4"|) shows a clear A^-type dependence and is proportional to the color factor Cf from the 
rescattering between an antiquark (or a quark) and a gluon. 



III. TRANSVERSE MOMENTUM BROADENING IN HEAVY QUARKONIUM PRODUCTION 

In this section we use the same technique reviewed in last section to calculate the transverse momentum broadening 
of heavy quarkonium production in both hadron-nucleus and nucleus-nucleus collisions. 

The heavy quarkonium's transverse momentum broadening in hadron-nucleus collision was often attributed to the 
initial-state multiple scattering between the active parton of the projectile and soft partons of the nuclear target before 
the hard collision to produce the heavy quark pair [jjjj . Calculation of such initial-state rescattering should be very 
similar to that for the Drell-Yan broadening, except that the quark-antiquark annihilation is accompanied by a much 
larger gluon-gluon fusion subprocess. If one considers only the gluon-gluon fusion subprocess, one should expect to 
have Eq. (|14|) for the heavy quarkonium broadening with the overall color factor C p replaced by C 'a — N c — 3 due to 
the difference in color factors between gluon rescattering and quark rescattering. The initial-state rescattering alone 
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leads to the naive expectation for the ratio of broadening between heavy quarkonium and Drell-Yan as Ca/Cf = 2.25, 
which is much smaller than the data [151 ]. 

However, as discussed in the Introduction of this paper, the net broadening of heavy quarkonium's transverse 
momentum in hadron-nucleus collision is a combined effect of the initial-state interaction and final-state rescattering 
between the produced heavy quark pair and the nuclear matter. This is because a heavy quarkonium is very unlikely 
to form at the same time when the heavy quark pair was produced. Since the final-state rescattering is sensitive 
to the detailed dynamics that transmutes a heavy quark pair into a bound quarkonium, we calculate the final-state 
contribution to heavy quarkonium broadening in both NRQCD and Color Evaporation models. 




(a) 



FIG. 3: Lowest order Feynman diagram for light quark-antiquark annihilation (a) and for gluon-gluon fusion to a pair of heavy 
quark. 



A. Color Evaporation Model 



In the Color Evaporation model (CEM), heavy quarkonium production is factorized into two steps: the production 
of a pair of heavy quarks with an invariant mass Q followed by a non-perturbative hadronization process with an 
universal transition probability for the pair to become a bound quarkonium 20]. It was assumed that the transition 
probability is the same for all heavy quark pairs whose invariant mass is less than the mass threshold of producing 
two open flavor heavy mesons, and the cross section for producing a heavy quarkonium, H , can be written as !2ij] 

OEM _ rp / Q ,r,2 aa hA^QQ 



«h = F QQ ^ H / dQ* > (15) 

where Fqq_^ h is a non-perturbative transition probability and is independent of the color and angular momentum 
of the heavy quark pair. There is one transition probability for each heavy quarkonium state, H . In Eq. (|T5j). the 
inclusive cross section for producing a pair of heavy quarks of invariant mass Q can be factorized as [35| 

= A^J Mta/HV) J dx<t> b/A { X ) , (16) 

a, b 

where ^ a b sum over all parton flavors, and da ab ^QQ j dQ 2 is a short-distance hard part for two partons of flavor 
a and b to produce a pair of heavy quarks of invariant mass Q. At the lowest order, they are given by the light 
quark-antiquark annihilation and gluon-gluon fusion subprocess, as sketched in Fig. [3] The transition probability in 
Eq. (|15p is assumed to be universal and independent of how the heavy quark pair was produced. It fixes the overall 
normalization for the cross section of heavy quarkonium production in different collision processes and provides the 
predictive power of the model. The model has been reasonably successful when comparing with data of inclusive 
heavy quarkonium production (32l [33j . 

Similar to Fig. [2j the leading order double scattering diagrams for producing a heavy quark pair are sketched in 
Fig. U for quark-antiquark annihilation subprocess, and in Fig. [5] for gluon-gluon fusion subprocess, respectively. The 
blob in the quark-antiquark annihilation subprocess in Figgis given by the diagram in Fig. [3]Ja), and the blob in the 
gluon-gluon fusion subprocess in Fig. [5] is given by a sum of the three diagrams in Fig. [3fb). 

In CEM, the transverse momentum broadening of a heavy quarkonium is equal to the transverse momentum 
broadening of the parent heavy quark pair, since the transition probability from a heavy quark pair to a bound 
quarkonium is given by a constant, Fqq^ h . We use the same method reviewed in last section to calculate the 





(a) 




FIG. 4: Leading order double scattering diagrams for qq 
scattering (b), (c), (d), and (e). 



QQ: initial-state double scattering (a), and final-state double 




FIG. 5: Leading order double scattering diagrams for 
scattering (b), (c), (d), and (e). 



(a) 







(d) " (e) 

QQ: initial-state double scattering (a), and final-state double 



transverse momentum broadening of the heavy quark pairs. Similar to Eq. ([7]) in the Drell-Yan case, we have 



dQ 2 dqf r 



d 1r It jt^jzt- = X! / dx '^q/hi x ') [ dxdx 1 dx 2 [T F ^(x,xi,X2,p)H^ Q ^(x,Xi,X2,p,q,x'p') 
q J J 

+ T F F q > ( x > x i, Xi 2^p) H< o^Lqq ( x , x i, x 2, P, q, x'p')} 



+Tf r J(x,x 1 ,x 2 ,p)H^ ' (x,x 1 ,x 2 ,p,q,x'p')\ , (17) 



gg^QQ' 

where the superscripts, "(/)" and "(F)", indicate the initial-state and final-state rescattering, respectively, and the 



matrix element T F ' is given in Eq. © . T F F is given by 



T F I p(x,Xi,X2,p) 



<h_ dVx_ t^2_ elXlP +yl i(x-xi)p+y- e -t(x-x 2 )p+y 2 

2ir 2tt 2tt 

x ( PA | F a + )F° + (0)F+ fo)F+ a (y-)\PA). 



(18) 



( F) ( F) 

The matrix elements with final-state rescattering, T F and T FF , have the same expressions as corresponding matrix 
elements with initial-state rescattering, since the field operators in the definition of the multiparton matrix elements 
in the collinear factorization approach commute on the light-cone (3 11 ] . 

The partonic parts, H i n ~, are given by the Feynman diagrams in Fig. [3] with the quark line from hadron (top) 

qq 

traced with (7 • p')/2, the quark line from nucleus (bottom) traced with (7 • p)/2, and gluon lines contracted with 

(/) ' 



p a pP. The diagram with initial-state rescattering in Fig. |3Ja) contributes to H q ^QQ as 



H 



(i) 

qq^QQ 



H 



jUx) _ 87r 2 q s 
qq^QQ ~ N 2 — 1 



2tt x\ — x — it x 2 — x + it 



do 



dQ 2 



(19) 



where the lowest order partonic cross section from qq annihilation to a heavy quark pair of invariant mass Q is given 
in Ref. [H]]. All four diagrams with the final-state rescattering, in Figs.|3Jb), (c), (d), and (e), contribute to H-\^qq 
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H (F) = H ^ 

qq^QQ qq^QQ 



8n 2 a s 



i 



C A 



2n xi 



dQ 2 



(20) 



The final-state contribution in Eq. (f2"0"|) is very similar to the initial-state contribution in Eq. (fit?)) except the overall 
color factor and the location of the unpinched poles. The difference in the location of the pinched poles, indicated 
by the sign difference of the ie, is a consequence of the order of the rescattering taken place either before or after the 
hard collision. The overall color factor for the final-state interaction, Ca in Eq. (|2H)) . indicates that as far as the color 
is concerned, the rescattering of a heavy quark pair is effectively the same as that of a color-octet gluon when the 
effect is calculated in the Color Evaporation model. 

Similarly, we evaluate the double scattering diagrams from gluon-gluon fusion subprocess in Fig. O and obtain their 
contribution to the partonic hard part, 



H 



(i) 

~gg- 



= H 



for the initial-state double scattering, and 



1 



C A 



1 



1 



1 



2ir xi 



x — ie X2 



ie 



do 



gg^QQ 
dQ 2 



H 



(F) 

gg^QQ 



rUMMM 

gg^QQ 



8n 2 a s 
N 2 -l 



C A 



2ir x\ 



ie x 2 



do 



gg^QQ 
dQ 2 



(21) 



(22) 



for the final-state double scattering. We find that the contribution to the gluon-gluon fusion subprocess from the 
final-state double scattering is the same as that from the initial-state interaction. That is because the rescattering 
effect of a heavy quark pair is the same as that of a color-octet gluon when the effect is calculated in the Color 
Evaporation model. 

Substituting the partonic hard parts in Eqs. (Q~9|), (120[) , (|2"Tj) . and (|2"2"|) into Eq. (fTT|) , integrating over the momentum 
fractions, x\ and x 2 of the rescattering gluons under the leading pole approximation, we obtain the leading double 
scattering contribution to the g^-moment of producing a heavy quark pair in hadron-nucleus collisions, 



do D 

,22 hA—>QQ 
dqT<lT dQ 2 dq 2 



87r 2 a s 
N*-l 



/ dx'^ /h (x') / dx C p tW a (x)+C a tW{x) 

q J J 

do 



do. 



+ I dx'<t> g/h {x') / dx C A T {I /A (x)+C A T ( Jl(x) 



qq^QQ 
dQ 2 



dQ 2 



(23) 



where the quark- gluon correlation function, T^ A , is given in Eq. the T^,' A is given by 



-,(-0 



dy 
2ir 



dyi dy 2 
2tt 



9(y- ~y-)9(-y-) 



x ZZT (pA\F a + (y^F^(0)F+(y-)F+ a (y^\p A ) 



xp~< 



(24) 



and T ( J ] A and T^fj 
are replaced as 



are given by the same expressions in Eq. (fTT|) and Eq. (|24|) . respectively, except the 9- functions 



9(y-- y -)d(-y 2 )^e(y^-y-)d(y^), 



(25) 



due to the different order of rescattering (or the sign of the ie for the unpinched poles). 

By integrating over the invariant mass of the heavy quark pair we can define the heavy quarkonium transverse 
momentum broadening in CEM as 



/HQ 



t 2 2 

dq T q T 



Ulf, 



47/17 



dQ 



do D 

2 a °hA~,QQ 

dQ 2 dq^ 



'4 m 7 



dQ 



2 d< ?hA^QQ 

dQ 2 



(26) 



where ttlq and Mq are heavy quark mass and open flavor heavy meson mass, respectively. 

As pointed out in Ref. (3lj . the field operators on the light-cone in the definition of the multiparton matrix elements, 
as those in Eqs. (jlip and (|24[) commute. The matrix element with initial-state rescattering is equal to corresponding 
matrix element with final-state rescattering, if the phase space interaction of these two matrix elements, such as the 
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J dy~ dy^[ dy^Oiy - — 2/i~)^( — y^) in Eq. (fTTj) for initial-state rescattering and J dy~ dy± dy^OiVi — for the 

corresponding final-state rescattering, are the same (l3| . However, the phase space integration for the final-state 
interaction in heavy quarkonium production may not cover the full size of the nuclear medium if the heavy quark 
pair becomes a physical quarkonium or transmutes into a color singlet pre-hadron quarkonium state before the pair 
exits the nuclear medium. Rescattering between a quarkonium and nuclear medium and that between a colored heavy 
quark pair and the same medium could be different, and lead to a different heavy quarkonium broadening. 

It was argued in Ref. [HI that a physical quarkonium state is likely to form outside nuclear matter in hadron- 
nucleus collision. Therefore, the matrix elements with final- and initial-state rescattering could have the same phase 
space interaction. If we assume that the matrix elements with the final- and initial-state rescattering are the same, 
qfl ( x ) = q/A ( x ~) ano ~ ^gfl ( JC ) = g'/A. ' ana ~ assumc the same model for both twist-4 quark-gluon and gluon-gluon 
correlation functions, 

Tff A (x)=\ 2 A*/ 3 ct> g/A (x), (27) 
we can express the heavy quarkonium broadening in hadron-nucleus collisions as, 

A(« M = ( |fe A 2 A^) (C F + C A )a q(l + 2C A a gg (2g) 

N„ 1 / (7 qq + Ugg 



where a qq and a gg are the lowest order inclusive cross sections from the qq — ► QQ — > H and gg — > QQ — ► H 
subprocess, respectively. They are given by [H, [34| 

r iM Q r r da - - 

*w = Fqq^h J a 2 dQ 2 J2j dx^ q/h (x') j dx<f> q/A (x) q J Q2 QQ (29) 

r iM Q r r da - 

° 9a = F QQ^H J 4 2 dQ 2 j dx'4> g/h (x') J dx<j> g/A {x) , (30) 

where the lowest order partonic cross sections to produce a pair of QQ are given in Ref. (34|. If the gluon-gluon 
subprocess dominates the heavy quarkonium production rate, a gg 3> a qq , we have 

A(4)g™«2^(^A 2 ^ . (31) 

By comparing the Drell-Yan broadening in Eq. (TJJI and the leading heavy quarkonium broadening in Eq. (|3"TT) . we 
conclude that the leading contribution to heavy quarkonium transverse momentum broadening in hadron-nucleus 
collision, calculated in CEM, is about 2C a /Cf = 4.5 times Drell-Yan broadening. 



B. Non-Relativistic QCD Model 



The application of NRQCD to the production of a heavy quarkonium H in hadronic collisions relies on the proposed 
factorization formalism I19II. 



NRQCD 
a hA^H 



f dx (j) a / h {x) J dx(j) b/A {x) 



ab->QQ[n] 



(O h (u)) 



(32) 



where a and b run over all parton flavors, and -ff a 6^QQ[n] are perturbatively calculable coefficient functions for produc- 
ing the heavy QQ[n] states. The state of the heavy quark pair, [n], is characterized by the pair's rotational, 2s+1 Lj, 
and color quantum numbers. The coefficient function for producing each QQ[n] state is perturbatively calculable in 
a power series of the strong coupling constant, a s . The matrix elements of (O h (n)) in Eq. (|3"2"]) describe the non- 
perturbative hadronization dynamics and give the probability for the pair to become a physical heavy quarkonium 
H [H, (36[. The matrix elements should be universal. That is necessary for the predictive power of the NRQCD 
formalism. The expansion in Eq. (I32p is organized according to the effective power of the heavy quark pair's relative 
velocity. Although it still lacks a fully compelling proof for the NRQCD factorization formula in Eq. (13"2")) f36| . the 
formalism for heavy quarkonium production has had many successes, in particular, its success in interpreting the 
CDF data on J ftp and if)' production as a function of transverse momentum 0, [37| • 
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In NRQCD model of heavy quarkonium production, the transition probability from a heavy quark pair to a bound 
quarkonium is sensitive to the pair's rotational and color quantum numbers. Partonic multiple scattering has a 
potential to change not only the heavy quark pair's momentum, but also the pair's color and other quantum numbers. 
Therefore, the transverse momentum broadening calculated in NRQCD model is not necessary the same as that 
calculated in CEM in last subsection. If the difference is significant, a precise measurement of transverse momentum 
broadening could shed some lights on heavy quarkonium's production mechanism. 

We assume that the g^-moment of heavy quarkonium production, J dq^(q 2 ^) n da^ < ^P /dq 2 ^ with n > 0, can 
be factorized in the same way as the th -moment given in Eq. (f32|) . We calculate the leading double scattering 
contribution to the broadening of the first moment. Like the calculation done in CEM in last subsection, the partonic 
double scattering diagrams are given in Figs. [4] and [5] Similar to Eq. (fl~7|) . we obtain 



/ 



j 2 2 aa hA^H 
dq T q T ^ 



E 



dx'4>q/ h (x) j dxdxx dx 2 [T^(x,xi,x 2 ,p)H^ H (x,xi,x 2 ,p,q,x'p') 

+ T F F q } {x,xi,x 2 , p)H^ H (x,xt,x 2 ,p,q, x'p')] 
dx'(f> g / h (x') J dxdxxdx 2 [T^l(x,Xi,x 2 ,p)H g P> _^ H (x,xi,x 2 ,p,q,x'p') 



+ t ff ( x > x i, x 2, P)H ( g g^ H {x,Xi, x 2 ,p, q, x'p')] , 



with the partonic cross sections defined as 

H$£> H (x,x u x 2 ,p,q, x'p') = £ H% F) (n) (O h (n)> , 



(33) 



H^h (x,x 1 ,x 2 ,p,q, x'p') = £ Hi 1 /) („) { o H (n)) , 



(34) 



where ^ n sums over all possible QQ states, n, with appropriate spin and color quantum numbers [38j |. and H^' F \n) 



and Hgg (n), whose dependence on parton momentum fractions and kinematic variables are suppressed, repre- 
sent partonic hard parts for producing a heavy quark pair at a quantum state n from quark-antiquark annihilation 
subprocess and gluon-gluon fusion subprocess, respectively. 

The partonic parts for the quark-antiquark annihilation subprocess, Hgg (n) in Eq. (I34|) . are derived from diagrams 
in Fig. 21 The single diagram in Fig.[4ja) contributes to Hgl\n). Under the leading pole approximation, the initial- 
state rescattering does not change the nature of the s-channel qq — » QQ subprocess, which produces a heavy QQ pair 



in a color octet and spin-1 state: n — 6 s\ . The corresponding hard part is given by 



8n 2 a s 
N?-l 



1 



1 



1 



2tt x-i 



it x 2 — x 



where ( 3 s[ 8 ^) is the lowest order short-distance coefficient for qq 



tt(0)/3c(8)\ 
n qq \ °1 I i 



QQ( 3 S^) subprocess and is given by 



n 3 a 2 16 



(35) 



(36) 



with the mass of a quarkonium: M = 2mg. The final-state rescattering effect to the qq annihilation subprocess 
comes from the four diagrams in Figs. HJb), (c), (d), and (e). The additional gluon rescattering in the final-state 
allows more quantum states for the produced QQ pair. At this order, we have the following nonvanishing states: 



3o(8) 3p(l,8) 



i q\°> 3 p 



j =1 2 , and corresponding hard parts, 

8ir 2 a. 



ff(f)/3p(l,8) s 
n qq \ r J=\,l) 



N 2 -l 
~~ N 2 -l 



C A 



1 



1 



2n x\ — x - 
1 



it x 2 



le 



2ir x\ — x + it x 2 — x — it 



tt(1)/3 p (1,8) s 
n qq \ r .]=l,2) > 



(37) 



where H$( S S[ 8) ) is given in Eq. (|36[) and the leading order P-wave contribution is given by 



<( 3 if } ) 



5 W 
3 H ™ 



(D ( 3p 2 (8) 



S(s-M 2 ). 



TT A aj 40 1 
M 3 27 3m 2 Q 

ir 3 a 2 64 1 c/ , 



(38) 
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The nonvanishing contribution to the 3 Pj]l^ > 2 s t a t es is a consequence of the gluon rescattering, which effectively 
provides a gg — > QQ subprocess. 

Similarly, partonic parts for the gluon-gluon fusion subprocess, Hgg (n) in Eq. (|34p . are derived from diagrams in 
Fig. [5l Unlike the quark-antiquark annihilation subprocess, the heavy quark pair produced in gluon-gluon subprocess 
can have more than one quantum state. For the initial-state rescattering diagram in Fig.[5ja), the heavy quark pair 



with n = 1 Sq 1 ' S \ 3 Pj 1 = ^2 can au §i ye nonvanishing contribution to Hgg 1 (n). The four final-state rescattering diagrams 



(I) 



in Figs.^b), (c), (d), and (e) can produce heavy quark pairs with n 
from gluon-gluon fusion diagrams in Fig. [5] as 



,(8) 



and 3 P 



(8) 



,7=0,2- 



We obtain the hard parts 



rr(/)/lo(l,8) 3 p (l,8) \ 
^gg \ °0 ! 0=0,27 



rr(F) ( lo(8) 3p(l 
n qq \ D ' C 



7=0,27 



N?-l 



N2-1 



C A 



C a 



1 


1 




1 


2-7T X\ 


— x — 


ie X2 


— x + ie 


1 


1 




1 


27T X\ 


-x + 


ie x 2 


— x — ie 



F/-(0Vle(l>8) 3n(l,8) n 

gg \ °o > 0=0,2/* ' 

H (P)(lcW 3p(8) n 
39 V J i 0=0,2/ ' 



(39) 
(40) 



where Hq^^Sq 1 ' , 3 -Pj=q^) are trie lowest order cross sections without the rescattering, 

i 

-B„J(5 - M 2 ) , 



rr(0)/lc( 
-"<J9 I ^0 



4 m ) 



M 3 
M 6 tuq 



M 3 



■B r , 



hm 2 Q 



S(s-M 2 ). 



(41) 



with a color factor B n = 2/9 for a color-singlet state and 5/12 for a color-octet state, respectively. 

Having obtained the short-distance hard parts, we can derive the leading double scattering contribution to the 
g^-moment of heavy quarkonium production in Eq. (|33[) . In this paper, we limit ourselves to the direct production of 
spin-1, 5-wave heavy quarkonia, such as J/i/j, T and etc. After neglecting the color singlet QQ states not equal to 3 Si, 
we are left with the following QQ states, which contribute to the direct production of a H( Si) heavy quarkonium, 



3c(8) 3p(8) \ 
3=1,21 



•i'i - wr.\ ■"./',- i,j • ^H( 3 Si), 

From Eq. ([34]), we obtain by summing over all QQ states in Eq. (|42]) . 



(42) 



H 



(i) 

qq^H 



Nl-l 
8ir 2 a s 



1 



1 



1 



2tt x\ 
1 



x - 
1 



ie X2 — x 
1 



c F H^es[ a) )(o H es[ B) )) 



2ir x\ 



ie X2 — x 



^F Oqq , 



Nl-1 

where the lowest order quark-antiquark annihilation cross section is defined as 

7r 3 a 2 16 



M 3 27 



S(s~M 2 )(O h ( 3 s[ s) )) 



(43) 



(44) 



with the nonperturbative NRQCD matrix element, (O h ( 3 s[ 8 ^)}, for a QQ[ 3 S'^ 8 ' ) ] pair to become a heavy quarkonium 
H. Similarly, we have the final-state rescattering contribution, 



H 



(F) 

qq^H 



8n 2 a s 
JV?-1 



87ra s 



1 



1 



2tt x\ — x 4- ie X2 



1 



1 



1 



(c A H${*si a) )(0 H esW)) 
+ff«( 3 pf)<0^ 

c A b 



2ir x\ — x + ie X2 — x — ie 
where we defined the P-wave contribution as 

ee ffW( 3 pf))<0^( 3 pf))) + J ffW( 3 pf )(0^( 3 Pf )) 
TT 3 a 2 80 



(45) 



M 3 27 



6(3- M 2 )(O H ( 3 P (8) )) . 



(46) 
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In deriving the last equation, we used the heavy quark spin symmetry |38j . 

(0 H ( 3 P ( f ] )) = (2J + 1)<0 H ( 3 P O (8) )) . (47) 

From Eqs. (|39[) and (|40[) , after we neglect the color singlet channels that have wrong spin and angular momentum 
quantum numbers, the gluonic contribution from the initial-state and final-state rescattering to the direct production 
of H( 3 Si) are effectively the same except the sign of the ie for the unpinched poles, 



H 



(i) 

gg^H 



H 



(F) 

gg^H - 



87r 2 a s 


' 1 


1 




1 


Ni-1 
8-K 2 a s 


2tt x\ 
' 1 


— X 

1 


— ie x 2 


— x + ie 
1 


Ni-1 


2n x\ 


— X 


+ ie x 2 


— x — ie 



Cao 



*(0) 



where the gluon-gluon fusion cross section is defined as 

(o H Csi s) )) + 4-<o*(»p < 8 >)) 



7 (0) 
99 



ir 3 a 2 5 
M 3 12 



8{s-M 2 ) 



(48) 



(49) 



In deriving Eq. (I48|) , the heavy quark spin symmetry was used to relate the matrix elements of different P-wave states 
to the matrix element of (O h ( 3 Pq 8 -*)). 

Substituting the partonic cross sections in Eqs. ([4*3"]) , l[35]). and (|4*5)l into Eq. (|3"3")) . and integrating over gluon 
momentum fractions, x\ and x%, by taking the residue, we obtain the double scattering contribution to the g^-moment 
of heavy quarkonium production in NRQCD model as, 



drr D 
dqTQT dq% 



8-K 2 a s 

W^i 



J2 ! dx'4> q/h {x') f dx \T§ A {x)C F a q ^ +T { J ] A {x) (c A a^+a 



+ / dx'<p g/h (x') / dx \T^ A (x)C A a^+T { J}(x)C A a^ 



(50) 



which has a very similar expression as that in Eq. (|2"3")) derived in CEM in last subsection. 

If we use the same model for the quark-gluon and gluon-gluon correlation functions as that used in CEM calculation 
in last subsection, we obtain the heavy quarkonium broadening in NRQCD model as 



A((? 2 ) NRQCD 



K 2 a s ^ 2 ^1/3^ ( C F + C A )cr q ^> + 2C A Ugg +a qq 



(i) 



r (0) (0) 

r qq ' "gg 



mu ; \N*-1 

where the leading order cross sections calculated in NRQCD model are given by 



(51) 



u qq 



^2 / dx ' fig/hi^') / dxcj) q/A {x)a q Y\ 

q J J 

a g °J = [ dx'^ glh [x') [ dx^ a/A {x)af g \ 



(52) 



with the partonic cross sections given in Eqs. (|44[> . (|46j) . and (|49p . respectively. 

From the transverse momentum broadening in Eq. (|2"8|) calculated in CEM in last subsection and that in Eq. (|5"Tj) 
calculated in NRQCD model, it is clear that the leading double scattering contribution to the broadening calculated 
in these two models have the same expression if one neglects the P-wave contribution in NRQCD approach. Since 
the P-wave contribution is smaller than the S'-wave contribution, and the gluon-gluon fusion subprocess dominates 
the heavy quarkonium cross section, we expect that A(<7f,)^Q QCD w A((7y)§Q M « (2 CU/Cf)A(<7^)dy- 

In both CEM and NRQCD approach to the production of quarkonia, H( 3 Si), we can also calculate the broadening 
effect on those quarkonia that were produced from the decay of either excited or high spin states of produced quarkonia, 
known as the feeddown mechanism of the quarkonium production. Since the g^-moment in Eq. |T]) is normalized by 
the cross section (the th -moment), and the rescattering takes place at the parton- level, we expect that the feed-down 
mechanism is not very sensitive to the quarkonium broadening while it is much more sensitive to the quarkonium 
production rate. We will come back to the role of the feeddown mechanism in quarkonium broadening in Sec. IVlwhen 
we present our numerical results. 
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C. Transverse momentum broadening in nucleus- nucleus collisions 



In this subsection, wc extend our calculations of heavy quarkonium's transverse momentum broadening in hadron- 
nucleus collisions to the broadening in nucleus-nucleus collisions. We discuss the similarities and differences between 
the hadron-nucleus and nucleus-nucleus collisions, and the role of transverse momentum broadening in probing the 
properties of the dense and hot QCD matter created in high energy nucleus-nucleus collisions. 



A 



B 



(a) 



(b) 



FIG. 6: Sketch of heavy quarkonium production in hadron-nucleus collisions as viewed in the target rest frame (a), and that in 
nucleus-nucleus collisions as viewed in the center-of-mass frame (b) . The thin and thick lines indicate the incoming parton and 
the outgoing heavy quark pair, respectively. The cross indicates potential rescattering point with soft partons of the nuclear 
medium. 



The broadening of heavy quarkonium's transverse momentum is a consequence of both initial-state and final-state 
rescattering in nuclear medium. In high energy hadron-nucleus collisions, both initial-state and final-state rescattering 
probe the same properties of a normal nuclear matter. If the scattering process is viewed in the rest frame of the 
nucleus, as sketched in Fig. [H^a), the incoming active parton and the produced heavy quark pair move very fast along 
the direction of the incoming hadron and interact with only partons of the nucleus near the same impact parameter. 
This picture supports our approximation that the correlation functions for final-state interaction are about the same 
as that for initial-state interaction. If the same process is viewed in the center of mass frame, we get the same 
conclusion that only nuclear partons near the same impact parameter can participate in the rescattering because the 
nuclear matter is moving very fast in this frame. 

However, in high energy nucleus-nucleus collisions, the situation can be very different. If, other than the hard 
parton-parton scattering to produce the heavy quark pair and corresponding rescattering discussed above, other 
partons from two colliding nuclei did not interact in the collision, the transverse momentum broadening of heavy 
quarkonium production would be a simple superposition of the broadening in hadron-nucleus collision. In this picture, 
the leading contribution to the broadening in a collision between a nucleus A and a nucleus B would be given by the 
same expression in Eq. (pJTj) with the A 1 / 3, replaced by A 1 / 3 + B 1 ! 3 or more precisely by Lab/i~o, where ro ~ 0.8 fm 
and Lab is an effective medium length in nucleus-nucleus collision (39j | . The value of ro is obtained by letting 
L p a ~ tqA 1 / 3 w (3/4:)Ra with the nuclear radius Ra- 

But, as indicated by the early data from RHIC, soft partons from two nuclei do interact to form a dense and 
hot QCD quark-gluon medium in high energy nucleus- nucleus collisions. The final-state rescattering between the 
produced heavy quark pair of transverse momentum and the almost stationary or slowly expanding hot QCD 
medium in the center of mass frame of nucleus-nucleus collision, as sketched in Fig.[S£b), is very unlikely to broaden 
the pair's transverse momentum. Instead, the final-state interaction could suppress the production rate of the leading 
(or large momentum) colored and coherent heavy quark pair due to the energy loss (23j . which is responsible for 
the observed jet quenching (or the suppression of leading hadrons or heavy quarks) at RHIC T], and decrease the 
averaged transverse momentum (q^). On the other hand, the initial-state interaction in nucleus-nucleus collisions is 
likely to broaden the transverse momentum of the active parton in the same way as that in hadron-nucleus collisions. 
Therefore, the measured (qf.) in nucleus- nucleus collision is a consequence of two competing effects: the initial-state 
interaction tries to broaden the transverse momentum while the final-state rescattering in a slowly expanding medium 
tends to reduce the pair's transverse momentum. The detailed analysis and calculation of the competing final-state 
effects on the quarkonium broadening in nucleus- nucleus collisions requires a careful modeling of the hot and dense 
medium and is beyond the scope of this paper. 

Precise measurements of transverse momentum broadening of heavy quarkonium production in relativistic heavy 
ion collisions should provide very valuable information on the formation of the dense and hot quark-gluon medium 
and its properties. In nucleus- nucleus collisions, a deviation of the transverse momentum broadening from the simple 
superposition of that measured in hadron-nucleus collisions clearly indicates a change of nuclear matter properties 
from the interaction between soft partons of colliding nuclei. It should indicate the formation of a dense quark-gluon 
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medium when the measured transverse momentum broadening is equal or less than the expected broadening from the 
initial-state interaction alone. 

If the long-range interaction of soft partons from two colliding nuclei is quick and strong, and the dense quark-gluon 
medium is formed very early in relativistic heavy ion collisions, the initial-state interaction in nucleus- nucleus collisions 
could be different from a superposition of the initial-state effect in hadron-nucleus collisions due to the modification 
of nuclear matter. In order to independently test the initial-state effect from the final-state rescattering, we calculate 
the transverse momentum broadening of Z (as well as W) bosons in relativistic heavy ion collisions at the LHC in 
next section. 



IV. TRANSVERSE MOMENTUM BROADENING OF Z (AND W) PRODUCTION AT THE LHC 

The lack of the final-state interaction of a Z (or a W) boson when it is reconstructed from its leptonic decay 
channels makes its transverse momentum broadening in high energy nuclear collisions an ideal probe for the initial- 
state interaction, in particular, the density of nuclear medium in the early stage of relativistic heavy ion collisions [25j . 
If the long-range soft gluon interactions between two incoming heavy ions were so strong and a dense nuclear medium 
was formed before the short-distance creation of a Z (or a W) boson, the transverse momentum broadening would 
be a clean probe of the density of the dense medium. Otherwise, the transverse momentum broadening of Z (or W) 
bosons in nucleus-nucleus collisions would be a simple superposition of that in hadron-nucleus collision. Therefore, by 
measuring the broadening of Z (or W) bosons in both hadron-nucleus and nucleus-nucleus collisions, we could learn 
valuable information on whether the dense quark-gluon medium could be formed at a very early stage in relativistic 
heavy ion collisions ■ 




(b) 



FIG. 7: Leading order double scattering diagrams for production of a Z (or a W) boson in hadron-nucleus collisions via the 
initial-state interaction (a) and the possible final-state rescattering if the vector boson is reconstructed from its hadronic decay. 

In Fig. [3 we sketch the leading double scattering diagrams for the Z (or W) production in hadron-nucleus collisions. 
The diagram in Fig. [7ja), which is almost identical to that for the Drell-Yan transverse momentum broadening, 
represents the initial-state interaction, while that in Fig. [T^b) represents the possible final-state rescattering if the 
vector boson could be reconstructed from its hadronic decay channels. For the Z and W bosons reconstructed from 
their leptonic decay channels, their transverse momentum broadening is mainly caused by the initial-state interaction, 
just like the broadening of the virtual photon in the Drell-Yan production. From the leading double scattering diagram 
in Fig.(7Ja) and following the same derivation for the Drell-Yan broadening in Sec. [ill we obtain the leading transverse 
momentum broadening of Z (or W) bosons in hadron-nucleus collisions as, 



= c f 



N*-l 
8n 2 a s (M w ) 



i 



\ 2 Al/3 



(53) 



where A| (A^) is the nonperturbative scale for the double scattering production of Z (or W) bosons. It is defined in 



the same was as that in Eq. (|13p and represents a ratio of nuclear four parton correlation function over a corresponding 
normal parton distribution. It is proportional to an averaged gluon field strength square in nuclear matter, (F +a F + a ) 
[40| . The Xz should be proportional to the typical transverse momentum or effective virtuality of soft gluons partici- 
pating in the double scattering. As we will discuss in next section, the value of the X\ (or Xyy) could depend on the 
momentum exchange of the hard collision, Q ~ Mz (or Mw), as well as the collision energy. 

If we assume that the dense quark-gluon medium in relativistic heavy ion collisions is not formed before the creation 
of the heavy vector boson Z (or W), we could apply our formula in Eq. (j53"j) to the broadening in nucleus-nucleus 
collisions as a simple superposition of the hadron-nucleus collision, and obtain the transverse momentum broadening 
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in nucleus- nucleus collision as 

A(^>^-^^^lA^, (54) 

where V = Z,W for the Z and W production, respectively. If the long-range soft gluon interactions between two 
colliding heavy ions are so strong that the dense quark-gluon medium was formed before the short-distance creation of 
a Z (or a W) boson, we expect the same formula in Eq. (fB"4"]) to be valid for the leading contribution to the transverse 
momentum broadening, but, with a new \ v different from that of a normal nuclear matter. The value of the effective 
Ay and its dependence on the collision geometry and collision energy should provide valuable informations on the 
formation of the QCD medium and its properties. 



V. NUMERICAL RESULTS 



In this section, we provide a numerical comparison between our calculated heavy quarkonium broadening in nuclear 
collisions with existing data from Fermilab and RHIC experiments, and make predictions for the transverse momentum 
broadening at the LHC. 

We have calculated in this paper the transverse momentum broadening of vector boson production in nuclear 
collisions in terms of the QCD factorization approach. We factorized the rescattering contribution to the broaden- 
ing into two parts: (1) the non-perturbative, but, well-defined universal parton-parton correlation functions, which 
represent the probability to find the scattering centers in the nuclear medium, and (2) corresponding parton-level 
rescattering subprocess, which are infrared safe and perturbatively calculable. As a result of the factorization, the 
normalization for the transverse momentum broadening from partonic double scattering is uniquely fixed by the size 
of the non-perturbative quark-gluon and gluon-gluon correlation functions. If we use the model in Eqs. (fl3|) and ([27]) 
to parameterize the correlation functions, the numerical results of our calculated transverse momentum broadening 
should be directly proportional to the value of the A 2 . 

The value of the A 2 , or more precisely, the value of the parton-parton correlation functions should not depend on 
which vector boson was produced. This is because of the universality and the long-distance nature of the correlation 
functions. However, the value of the correlation functions or the A 2 should depend on the momentum scale at which 
the hard part of the partonic scattering was evaluated. 

As explicitly demonstrated in last three sections, all phase space integrations for the partonic rescattering can be 
deformed into perturbative region, and both initial- and final-state parton-level rescattering are evaluated at a hard 
scale Q ~ 2toq or Mz (or Mw) for heavy quarkonium or Z (or W) production, respectively. From the model in 
Eq. (fT3")) (or (|27p). the non-perturbative parameter, A 2 , represents a ratio of nuclear four parton correlation function 
over a normal parton distribution. As demonstrated in Ref. by approximately decomposing a nuclear state into 
a product of nucleon states, the ratio (so as the A 2 ) can be expressed in terms of an averaged gluon field strength 
square, (F +a F^) . In this picture, the A represents the virtuality or the typical transverse momentum of the partons 
participating in the partonic rescattering, and certainly depends on the momentum exchange of the hard collision, Q, 
as well as the collision energy, y/s, which determines the available phase space for the collision. For hadron-nucleus 
collisions with a large momentum transfer, Q, we expect the A 2 to be proportional to ln(Q 2 ) from the parton shower 
leading to the hard collision. If the collision energy ^fs is very large and the momentum fraction of the active parton, 
x, is very small, we would expect the A 2 to be proportional to the saturation scale Q 2 S oc \/x s with S « 0.3 fill. 1421] . 
For the vector boson production, the typical momentum fraction of an active parton, x ~ Q/tJs. Therefore, we expect 
A 2 (Q)cxln(Q 2 )(y5/Q) 5 . 

We use data on the Drell-Yan broadening in hadron-nucleus collisions, which were measured by Fermilab experiments 
at the fixed target energy, y/s = 38.8 GeV [ll|, to extract the A^y From the value of Aq Y , we estimate the value of 
A 2 for producing a vector boson of invariant mass Q at a collision energy y/s as follows, 

ln(Q 2 ) / yfs/Q 



X (0) M Ady HiQ^)) U.8/(Qdy>J ' (55) 

with the yfs in unit of GeV, (Qdy) ~ 6 GeV, and 5 ~ 0.3. 

Fermilab experiments: E772, E789, and E866 have measured the transverse momentum broadening of the Drell- 
Yan, as well as J/ip, ip', and T production in hadron-nucleus collisions [HI, [22j. In Fig. [8j we plot the data on both 
the Drell-Yan and heavy quarkonium broadening as a function of atomic weight of nuclear targets. The broadening 
for the data was defined as a difference between the (^-moment in proton- nucleus and proton-deuteron collisions: 
A(<7^) = {q^) P A — (<1t)pD- By fitting the data on the Drell-Yan broadening as a function of A 1 / 3 — 2 1 / 3 , we obtain 
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FIG. 8: Data on transverse momentum broadening of heavy quarkonium as well as Drell-Yan production in hadron-nucleus 
collisions. Also plotted are theoretical calculations using Eq. (|28|l (solid lines) and Eq. (|51|l (dashed lines), derived by using 
CEM and NRQCD model, respectively. Three solid lines (from the top to bottom) correspond to J/ip, T, and Drell-Yan, 
while three dashed lines represent ip', J /ip, and T production from NRQCD model. The quarkonium broadening calculated in 
NRQCD model is evaluated with (a) and without (b) quarkonia from the feeddown mechanism. 



Xjjy ~ 001 GeV 2 , whic h g ives the bottom solid line for the Drell-Yan broadening in Fig. and is consistent with the 
value extracted in Ref. [13| . 

In Fig. [HJa), we plot our theoretical calculations of transverse momentum broadening of direct heavy quarkonium 
production in hadron-nucleus collisions at the Fermilab fixed target energy, y/s — 38.8 GeV. To obtain the numer- 
ical results of theoretical calculations, we use CTEQ6L for nucleon parton distribution functions [H| and EKS98 
parametrization for nuclear parton distributions (nPDF) [44| to evaluate the leading order production cross sections 
in Eqs. and (fSTj) . The non-perturbative parameter, A 2 , in Eqs. ([28")) and (foTj) could be slightly different for 
J/ip and T production due to the difference in the scale of hard collision, Q ~ 2mg. Using Aj^y ~ 0.01 GeV 2 , 
M J/V , = 3.1 GeV, My = 9.5 GeV, we estimate from Eq. ^ that X] /ip « 0.008 GeV 2 and X% w 0.011 GeV 2 at 

y/s = 38.8 GeV. We use Eq. (fS"Tj) to evaluate the broadening for tp', J/ip, and T (the dashed lines) in NRQCD model. 
The NRQCD matrix elements are taken from Refs. [1, E^j. The small difference between the broadening of J/ip and T 
is caused by the relative size between the quark and gluon contribution due to different sizes of matrix elements and 
the range of parton momentum fractions. For the direct production, J/ip and ip' have almost the identical broadening 
as shown in Fig. [8ja). Similarly, we use Eq. ([28]) to evaluate the broadening for J/ip (top solid line) and T (middle 
solid line) in CEM with m Q = 1.5 GeV and Mq = 1.85 GeV for J/ip and m Q = 4.5 GeV and Mq = 5.3 GeV for T 
production, respectively. The transition probability, Fqq^ h , cancels between the numerator and the denominator in 
Eq. (|28p . The small difference between these two solid lines are again caused by the relative size of the quark and 
gluon contribution due to the slightly different ranges of parton momentum fractions for J/ip and T production. 

We also test the effect of transverse momentum broadening on heavy quarkonia produced by the feeddown mecha- 
nism. Using the partonic hard parts derived in Sec. IIIIB1 we calculate in NRQCD model the transverse momentum 
broadening of quarkonia in high spin or excited states, which then decay into the observed spin-1 and S- wave quarko- 
nia, H( 3 Hi). We plot the calculated broadening with this feeddown mechanism in Fig. [8]^b). Although the feeddown 
mechanism could provide as much as 40% of J/ip production rate, its net effect on the transverse momentum broad- 
ening is very small because of the fact that the broadening defined in Eq. ([1]) is normalized by the inclusive cross 
section. 

For a large nucleus with the atomic weight A > 10 2 , some of the produced heavy quark pairs could transmute to a 
color singlet pre-quarkonium state (or even a physical quarkonium) before exiting the nuclear matter. Comparing to 
a colored heavy quark pair, these color singlet states should have a weaker interaction with the nuclear matter and 
get less broadening in transverse momentum. Therefore, we expect the theoretical curves in Fig. [H] to be slightly less 
steep than what were shown when A > 10 2 . 
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From Fig. [8j we conclude that perturbative QCD calculations of the quarkonium broadening based on both CEM 
and NRQCD model give a good description of existing experimental data in hadron-nucleus collisions. The major 
difference between the heavy quarkonium and the Drell-Yan broadening is naturally explained by the role of final-state 
interactions. Although the production mechanism in CEM and NRQCD model is different, these two models of heavy 
quarkonium production predict almost the same functional form for the transverse momentum broadening, as shown 
in Eq. (f28|) and Eq. (j5"Tj) , respectively. Since the P-wave contribution is relatively small, as shown in Fig. [H these 
two models predict almost the same transverse momentum broadening. In addition, both models predict that J /if; 
and T have effectively the same broadening in hadron-nucleus collisions other than a small difference caused by the 
available phase space (i.e., the available range of parton momentum fractions). 
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FIG. 9: RHIC data on averaged transverse momentum square of J/ip production as a function of the number of participants 
[4(| . The top panel is for the J ftp's produced in the central rapidity region while the bottom panel is for those produced in 
more forward or backward region. Also plotted are theoretical calculations using Eq. (|56p . Solid lines for Au-Au collisions and 
dashed lines for Cu-Cu collisions, respectively. 



In Fig. [51 we plot PHENIX data on averaged transverse momentum square of J ftp production at RHIC energies as 
a function of the number of participants, N part p?j| . The top panel is for the 3 ftp produced in the central rapidity 
region with \y\ < 0.35, and the bottom is for those produced in the forward and backward region with 1.2 < \y\ < 2.2. 
We also plot our theoretical calculations of the transverse momentum square by using 
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where {q^)\ pp - exp is fixed by the data point from proton-proton collision in Fig. [5] and A(gy.)AB is given by our 
calculation. We evaluate A{q^) ab in nucleus- nucleus collisions with an effective medium length Lab as 
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(57) 



in CEM. Eq. ([ST]) is obtained from Eq. (|2"8"]l by replacing the A 1 / 3 by Lab/ Tq. We can calculate the broadening 
in NRQCD model by applying the same replacement to Eq. (|5"Tj) . At the RHIC energy, we obtain Aj/^(RHIC) 2 w 
0.013 GeV 2 from Eq. (|5"5"j) . We calculate the Lab in Glauber model and convert it to N part [39l ] . 

In nucleus-nucleus collisions, soft gluons from the colliding ions can interact even if the collision is not a head-on 
or a central collision. Such interaction in a non-central collision may not be strong enough to stop all soft partons 
to form an almost stationary or slowly expanding hot medium of quarks and gluons. It certainly can slow down 
some of the colliding soft partons to change the characteristics of the nuclear matter, which could alter the final-state 
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interactions. As a result, the final-state interaction between the produced heavy quark pairs and the modified nuclear 
matter in nucleus- nucleus collisions, as sketched in Fig. [Hfb), generates less transverse momentum broadening if it 
does not reduce the transverse momentum due to energy loss. In Fig. [9l we plot our calculations of J/ip transverse 
momentum broadening by using Eqs. (|56[) and (|57[) . and keeping only the contribution from initial-state rescattering. 
The solid lines are for the Au-Au collision, while the dashed lines arc for the Cu-Cu collision. Our calculations are 
consistent with the data in both rapidity regions. 

In the central Au-Au collision, a hot and dense medium is produced. As discussed above, the averaged transverse 
momentum could be reduced, instead of the broadening, due to the energy loss of the produced heavy quark pairs 
when they interact with the hot and slowly expanding medium. A more detailed study of the momentum shift of the 
heavy quark pairs in such a hot medium could provide a more accurate description of the data in central region, and 
help the extraction of medium properties. 

Transverse momentum broadening of Z (or W) bosons in high energy nuclear collisions could provide a clean mea- 
surement of initial-state interactions, and help isolating final-state rescattering effect in heavy quarkonium production. 
However, because of the heavy mass of Z (or W) boson, only the LHC has a chance to measure the broadening reliably 
0. 

We use Eq. ([55)) to estimate the A 2 for the heavy vector boson production at y/s = 5.5 TeV, the averaged nucleon- 
nucleon collision energy in relativistic heavy ion collisions at the LHC. We obtain A 2 (LHC) « 0.035, 0.05 and 0.05 GeV 2 
for the production of J/ip, T, and the heavy vector boson Z (or W), respectively. Although the T mass is much 
smaller than that of a Z boson, the T and Z have the same A 2 (LHC) for the transverse momentum broadening due 
to the larger available phase space for T production [47]. In Fig. [TUJ we present our predictions for the transverse 
momentum broadening of vector boson production at the LHC. Using the estimated X 2 Z ^ W (LRC) w 0.05 GeV 2 and 
Eq. (|53p. we evaluate the transverse momentum broadening of Z (and W) bosons reconstructed from their leptonic 
decays, and plot the predictions for hadron-nucleus collisions as a function of atomic weight of the nucleus, A, in 
Fig. UOf a). We also plot the expected transverse momentum broadening of J/ip and T production in hadron-nucleus 
collisions in Fig. fTOT a). The curves for heavy quarkonium broadening are evaluated by using A 2 y^,(LHC) w 0.035 GeV 2 
and A^ (LHC) w 0.05 GeV 2 , and Eq. ([25]) from CEM without contributions from the feeddown mechanism. Eq. (|5Tj) 
derived from NRQCD model gives the similar results. The heavy quarkonium broadening in Fig. UOf a) is much larger 
than that of Z (or W) bosons because of the additional final-state effect, and the difference in color factor and the 
strength of the strong coupling constant, a s {Q). 




FIG. 10: Transverse momentum broadening of Z and W (lower set of curves) in hadron-nucleus collisions (a) and nucleus- 
nucleus collisions (b) at y/s — 5.5 TeV as a function of atomic weight of nuclear beam and the number of participants, N part , 
respectively. Also plotted are predictions (upper set of curves) for the transverse momentum broadening of heavy quarkonium 
production in hadron-nucleus collision at the LHC energy (a); and the initial-state only broadening in nucleus- nucleus collision 
at the same energy (b). 

In Fig. llOf bh we plot the expected transverse momentum broadening of vector boson production in Pb-Pb collision 
at y/s = 5.5 TeV at the LHC as a function of the number of participants, N part 39]. We calculate the Lab in 
Glauber model with inelastic nucleon-nucleon cross section <J^ N = 70mb at the LHC energy and convert it to N part 
in Pb-Pb collisions (3{|. For the Z (and W) broadening, we use Az/w(LHC) 2 ~ 0.05 GeV 2 , the same value used 
for the broadening in hadron-nucleus collisions in Fig. UOf a). Since the transverse momentum broadening is directly 
proportional to the A 2 , which is proportional to the gluon strength in the medium, a deviation from the predicted 
curves in Fig. fTDT b) signals the formation of the hot quark-gluon medium before the creation of the heavy Z (or W) 
bosons. 
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For 3/ip and T production in Fig. Mb), we use A^(LHC) « 0.035 GeV 2 and A 2 - (LHC) « 0.05 GeV 2 , and 
Eq. (|28|) from CEM to evaluate the quarkonium broadening. Eq. (|5~lj) from NRQCD model gives similar predictions. 
The plotted curves for 3/ij) and T production in Fig. HOf b) include only initial-state interaction. As discussed in 
Sec. 1111 CI the final-state rescattering in nucleus-nucleus collisions is more likely to reduce the transverse momentum 
broadening due to the energy loss, if a slowly expanding medium was produced. Therefore, we expect the curves in 
Fig. rroT b) to represent the maximum broadening of 3/i/j and T that will be seen in relativistic heavy ion collisions at 
the LHC [24|. 

If we could reconstruct Z and W bosons from their hadronic decay (e.g., into two jets), which might be impossible to 
do in the LHC environment, the hadronic jets from Z and W decay will have to interact with the nuclear medium. The 
final-state multiple scattering could generate momentum imbalance between these two jets and effectively introduce 
an apparent mass shift for the Z and W boson [9j. Such mass shift would provide additional information on the 
properties of the hot quark-gluon medium. 

VI. SUMMARY 

In this paper we evaluated transverse momentum broadening of vector boson production in hadron-nucleus and 
nucleus-nucleus collisions in terms of parton-level multiple scattering. We argued that the broadening, A(qj>)AB = 
{<It) ab — {<^)hhi defined as a difference between the averaged transverse momentum square measured in nuclear 
collisions and that in hadron-hadron collisions, can be systematically calculated in perturbative QCD in terms of 
the collinear factorization approach !7j. We factorize the contribution to transverse momentum broadening into 
the calculable short-distance partonic rescattering multiplied by universal parton-parton correlation functions, which 
specify the probabilities to find coherent scattering centers in a nuclear medium. We explicitly show that the par- 
tonic rescattering diagrams are evaluated at a perturbative hard scale Q ~ 2tjiq. We verify the universality of the 
non-perturbative parton-parton correlation functions by fitting the data on both the Drell-Yan broadening and the 
broadening of J ftp and T production to clearly demonstrate the predictive power of the QCD factorization approach. 

For the Drell-Yan virtual photon and Z (or W) boson production, we evaluated transverse momentum broadening 
from purely initial-state multiple scattering. We discussed the scale dependence of the non-perturbative parameter, 
A 2 , and extrapolate it from its value at the fixed-target energy to the collider energies. We calculated the broadening 
of Z (and W) bosons in both hadron-nucleus and nucleus-nucleus collisions at the LHC. We evaluated the Z (and W) 
boson broadening in nucleus- nucleus collisions as a superpostion of its broadening in hadron-nucleus collisions, and 
argued that a deviation from our calculation is a clear signal indicating that strong interactions between soft partons 
of colliding heavy ions took place before the short-distance creation of Z bosons. That is, the transverse momentum 
broadening of Z (or W) bosons could be a clean and excellent probe of the early stage dynamics of relativistic heavy 
ion collisions. 

For 3/ip and T production, we demonstated that the final-state interaction between the produced heavy quark 
pair and the nuclear medium is very important in understanding the existing data. We calculated the broadening in 
both CEM and NRQCD model, and clearly showed that the two models generate a small difference in the broadening 
and the broadening has a weak dependence on the feeddown mechanism. That is, our results are insensitive to 
the details of the hadronization mechanism and perturbatively reliable. We found that the leading contribution to 
heavy quarkonium broadening in hadron-nucleus collisions is about 2Ca/Cf = 4.5 times the corresponding Drell-Yan 
broadening, which gives a good description of the existing Fermilab data. 

We argued that the role of the final-state interaction to the transverse momentum broadening in nucleus- nucleus 
collisions could be very different from that in hadron-nucleus collisions. In hadron-nucleus collisions, both initial-state 
and final-state interactions involve only soft partons of the nucleus near the same impact parameter, and therefore, 
provide similar contributions to the transverse momentum broadening as long as the heavy quark pair stays in a 
colored state. On the other hand, soft partons from two colliding nuclei could strongly interact to produce a slowly 
expanding quark-gluon medium. This new medium could interact with the produced heavy quark pair if the pair is 
in a color octet state, and the interaction could be very weak or vanish if the pair is in a singlet state. The final-state 
interaction with this new medium in nucleus-nucleus collision, as sketched in Fig. [5] could reduce (instead, to broaden) 
the pair's transverse momentum. 

We predicted the transverse momentum broadening of vector boson (3/i/j, T, and W/Z) production in hadron- 
nucleus collisions and a maximum broadening in nucleus-nucleus collisions at the LHC in Fig. [TTJ] A more detailed 
calculation of final-state interaction between the produced heavy quark pair and the dense and slowly expanding 
quark-gluon medium produced in nucleus-nucleus collisions could provide a better conncetion between the measured 
average of the transverse momentum square and the properties of the dense quark-gluon medium in relativistic heavy 
ion collisions [48j . Having a precise transverse momentum broadening could also shed some lights on nonperturbative 
dynamics that determines the transverse momentum distribution of vector boson production in low qt region [l(| • 
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Finally, we emphasize that our QCD factorization approach for calculating the heavy quarkonium broadening can 
be applied for studying the nuclear dependence of quarkonium cross sections in high energy nuclear collisions [5(, as 
well as the nuclear dependence of the quarkonium's rapidity and transverse momentum distributions if the transverse 
momentum qx is large enough 0, §]. The advantage of the QCD factorization formalism is that we may readily 
quantify corrections in powers of a s . 
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